Fluctuation spectrum of quasispherical membranes with force-dipole activity 
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The fluctuation spectrum of a quasi-spherical vesicle with active membrane proteins is calculated. 
The activity of the proteins is modeled as the proteins pushing on their surroundings giving rise 
to non-local force distributions. Both the contributions from the thermal fluctuations of the active 
protein densities and the temporal noise in the individual active force distributions of the proteins 
are taken into account. The noise in the individual force distributions is found to become significant 
at short wavelengths. 
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I. INTRODUCTION 

Biological membranes are mixtures of basically lipids 
and proteins that actively participate in many biologi- 
cal processes. The lipids, being the basic component, 
form a bilayer in which proteins are included. Many of 
these membrane proteins are molecular machines, which 
can use some readily available energy source, for instance 
light or ATP, to perform different tasks . A prominent 
example is the transportation of ions across a membrane 
against a chemical potential gradient. Due to this con- 
stant activity a biological membrane is a system out of 
thermal equilibrium. 

One way to observe non-equilibrium behavior of mem- 
branes is to study their shape fluctuations. Lipid mem- 
branes in their fluid phase are flexible structures undergo- 
ing thermally excited shape fluctuations that can be ob- 
served directly by video microscopy P|- Indirect studies 
of the fluctuations, using micromanipulation techniques 
to measure the amount of membrane excess area stored 
in the fluctuations, have already given evidence that acti- 
vated proteins incorporated in lipid membranes enhance 
these fluctuations 0, Q • Efforts are being put into ex- 
tending these studies to use video microscopy to obtain 
the fluctuation spectrum directly ||5|j. 

A theoretical model was proposed in , which explains 
the activity induced enhancement of the fluctuations as 
being a result of the active proteins pushing on their fluid 
surroundings. This pushing could arise when a protein 
is pumping material across the membrane, or when it 
is changing its configuration as part of an active cycle. 
This model will be called the force-dipole model here, 
because it was formulated mathematically as each pro- 
tein contributing with a force-dipole to the hydrodynamic 
equations of the fluids surrounding the membrane. The 
consequences of this model for excess area stored in the 
fluctuations of a planar membrane was calculated in 
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and from this calculation it was found that the model 
constitutes a possible explanation of the results of the 
micromanipulation experiments. 

A video microscopy study of quasi-spherical vesicles 
would test the force-dipole model further, because the 
fluctuations at different wavelength can be measured, 
and not just the sum over all wavelength giving the ex- 
cess area. However, it is not easy to calculate the pre- 
dictions of the model for the fluctuation spectrum of a 
quasi-spherical vesicle directly, because the equations of 
motion for the surrounding bulk fluids would then have 
to be solved in a spherical geometry including the pres- 
ence of the force-dipoles. 

In the force-dipole model was reformulated into an 
equivalent but calculationwise more manageable formu- 
lation, which was called the Gibbs formulation for the 
system. One purpose of the present paper is to show 
how this reformulation makes it possible to calculate the 
fluctuation spectrum of a quasi-spherical vesicle. 

A second purpose of this paper is to include the non- 
thermal noise of the activity in the calculation of the 
fluctuation spectrum. This was neglected in 6], where 
the magnitude of the force-dipole generated by each ac- 
tive protein was treated as being constant in time. How- 
ever, there does not seem to be good reasons to assume 
that the forces with which a protein is actively pushing 
on its surroundings is constant in time as the protein 
moves through one cycle of, say, pumping an ion across 
the membrane. Estimates performed for a planar mem- 
brane in , where the effect of active noise were treated 
independently of the activities time-averaged constant ef- 
fect on the dynamics, showed that the noise can have a 
significant effect. Here we will include both effects for 
force-dipoles, the time averaged and the corresponding 
temporal noise, and calculate the combined effect on the 
fiuctuation spectrum of a quasi-spherical vesicle. 

The presentation of the calculation will begin in sec- 
tionm where the near-equilibrium dynamics of the mem- 
brane is established. The dynamics is formulated as a set 
of Langevin equations in section IIIIL and the appropri- 
ate statistics of the thermal noise is worked out. This 
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is in turn based on a calculation of the equilibrium fluc- 
tuations performed in appendix ^ The activity is in- 
troduced in section IIVI where the Langevin equations 
for the membrane dynamics is modified to include the 
active force-dipoles. The main result of the paper, equa- 
tion (|57|l . is in section where the activity modified 
Langevin equations are solved for the shape fluctuations 
of the vesicle and the different contributions to the result 
are estimated and discussed. In section IVll it is briefly 
sketched what impact the inclusion of the non-thermal 
noise has on the analysis of the existing micropipette ex- 
periments. Finally, a conclusion is given in section rVIII 



II. NEAR-EQUILIBRIUM DYNAMICS 

In this section we will set up the near-equilibrium dy- 
namics following the formalism developed in [9|. The 
notational conventions regarding differential geometry is 
the same as in 0,13 • 

Our starting point is an excess Hclmholtz free energy 
F of the membrane which is 



F = dA [2kH^ + A (np+ - rip-) H 
Jm 



-a^{(b-l) + (j)V{np+/(b,np-/cb)] . (1) 

The first term represents the bending resistance of the 
membrane with k being the bending rigidity and H the 
mean curvature. The second term is a coupling between 
the mean curvature and the protein concentrations per 
membrane area, np+ and n^- , where -I- and — label the 
two possible orientations that an asymmetric transmem- 
brane protein in a membrane can have. The third term 
is a consequence of the observation that inhomogeneities 
arising from local compression of the area of the mem- 
brane relax much faster than bending deformations '10|. 
The term therefore has the form of a local constraint 
with being a Lagrange multiplier field The field 
(j), which is constrained to be 1, is a linear combination of 
the different density fields ha of molecules in the mem- 
brane, with each species (including both proteins and 
lipids) weighted by the area the individual molecules 
prefer to take up in the membrane 
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■ ^ a^UA 
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The unitless number a is a global quantity that mea- 
sures the mismatch between the area that the membrane 
prefers to have. 



^preferred ^ ^2 "'^ dA UA , 



and the actual area of the membrane 



A^ 



dA 



M 



(3) 
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It is explicitly 



-^preferred 

A 



(5) 



The average of over the area of the membrane is ther- 
modynamically conjugate to the area A, and thus this 
average can be understood as the tension of the mem- 
brane. The fourth term in Equation ^ is an effective 
interaction potential between the protein fields np± . The 
functional form is left unspecified, but it should include 
at least the mixing entropy for the protein fields np± , and 
could also contain terms related to other interactions be- 
tween the protein molecules. The 0-field in this term is 
introduced for convenience, because then the potential V 
will not be present in the force balance equation of the 
membrane, see Equations (O to © later. 

Having specified the free energy in Equation we 
can now follow the formalism of 9] and write down the 
dynamics of the membrane. At low Reynolds number the 
equation of motion for the shape field R is simply the 
condition of force balance of the membrane. Following 
@ we ignore internal dissipation in the membrane and 
write the force balance as 
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The force /^.^ is the elastic restoring force per area of the 
membrane 
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/,,,„ ^2Ha^ - AkH [H^ - K) 



A. 



K{2H^ ~K) riA- TrVj^nA , 



/rs,a =9aCr0 - AuAdaH , 



(7) 

(8) 
(9) 



where we have introduced the density difference field 
nA = np+ — Up- . This field has a natural partner in 
the field = f^p+ + rip- , which will be used later. The 
force densities are the stresses of the bulk fiuids on 
the membrane surface. At low Reynolds number the bulk 
fluids are governed by the Navier-Stokes equation 



r]V^v± - Vp=^ = , 



(10) 



together with the incompressibility condition V • v± — 
0. Here v± — v±{r,t) are the velocities of the bulk 
fluids with r labeling the position and t time, are the 
pressure fields and r/ the viscosity. Given a solution of the 
bulk hydrodynamics we can calculate the stress tensors 
of the bulk fluids as 



T=^ = -p=^l + 77 Vv± + {Vv 



(11) 



where I is the identity tensor. Then the forces T are 



±n- 7=" 



\r=R 



(12) 
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To actually find a solution of the bulk hydrodynamics we 
need boundary conditions for the bulk velocity fields v± 
at the membrane surface. The boundary condition that 
is usually imposed when dealing with viscous fluids is 
the so-called no-slip boundary condition, where the bulk 
fluid at the boundary is constrained to have the veloc- 
ity of the boundary material. In it was argued that 
for a multicomponent membrane the general form of this 
boundary condition was a linear combination of the flows 
of the different species of molecules in the membrane. If 
we let be the current for the density ua entering the 
conservation law [9l| 



DtUA + D^jI = , 
then we can write such a linear combination as 



'"±\r=R 



dtR 



A 



(13) 



(14) 



where the phenomenological constants are areas of 
molecular sizes. They should satisfy the constraint 
^±^^ — 1 such that it is possible to have uniform 
motion of the whole system. In it was argued that 
the precise values of the L±. can influence the distribu- 
tion of proteins in the membrane if the membrane is sub- 
jected to a bulk shear flow with velocities of the order 
kBT/{r]L^). However, the velocities in the bulk fluids 
induced by thermal noise are of the order ksT / [t])?), 
where A is the wavelength of the undulation in question. 
Thus when we are interested in undulations much longer 
than the molecular length scales we can ignore shear flow 
effects arising due to thermal fluctuations. This means 
that the detailed values of the are unimportant for 
the coupling of the bulk fluid dynamics to the dynamics 
of the protein densities in the present problem. Thus we 
are allowed to make a convenient choice that will simplify 
calculations. We will choose: — a"^/(l + a) such that 



(15) 



A reason that the boundary condition in Equation l|15|) 
is convenient can be found by multiplying H13|l by a'^, 
summing over A and using the constraint ip — \ to flnd 
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(16) 



where w^ta — — dtR. Thus the motion of the mem- 
brane shape R is related to the motion of the bulk fluids 
v± at the membrane surface, without the protein densi- 
ties ripi entering directly. 

From Equations to © we see that the dynamics 
of the shape field R is coupled to the density difference 
field ua • The equation of motion for ua is a conservation 
law, Equation H13|l. but we need the expression for the 
current j^. This current can be divided into two contri- 
butions, — riAw" + Ja d' where the first term is the 



convective flow with v" = '^A^^jA/Cl2B''^^''^B) rep- 
resenting the kinematic velocity of the flow in the mem- 
brane and being the mass of a molecule of species A. 
The second term is the dissipative part of the current. 
Close to equilibrium we can write down a constitutive 
relation coupling this current linearly to possible driving 
forces. The only driving force we have that transform 
in the same way as under symmetry transformations 
(including the approximate symmetry of exchanging np+ 
and rip-) is the gradient of the corresponding chemical 
potential 



A 



SF 



Sua 



dV 



R,n^ ,0—1 



duA 



AH 



Therefore we can write the constitutive relation as 



A 



(17) 



(18) 



where JIaa is a phenomenological parameter related to 
the diffusion constant of the proteins, see Equation H65() 
later. 

Having worked so far with equations that are true for 
any membrane shape, we will now restrict ourselves to 
a shape which is fluctuating around a spherical aver- 
age shape with radius Rq, and only keep terms up to 
flrst order in deviations from this average shape. We ex- 
pand in spherical harmonics yim = 3^£m(^, 0), which are 
eigenfunctions of the Laplace operator on the sphere with 
eigenvalues —£{£+ 1)/Rq, and they are orthonormal 
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sin6'd6' yg,^,yem = StiSm', 
The expansion will be written as 

R = Ro + noRo ^ U^m3^£m , 

£,771 

riA — ^ nA,i7nyi77i , 

llT. = no,E + ^ n^,lmytm , 
£,rn 

(^0,<p + ^ O-^^lmyi-, 



(19) 

(20) 
(21) 
(22) 
(23) 



£.7n 



where the subscript indicates that the quantity has the 
value of the average state and ugm, nA,em, n^,im and 
o'ip/m represent small deviations from this state. Note 
that we follow ^ in making the simplifying assumption 
that the proteins are distributed symmetrically between 
the two possible orientations in the membrane such that 
no, A vanishes. 

The solution of the equations of motion to zeroth or- 
der in the deviations from a spherical membrane shape is 
that Ro and no,s are constant due to volume and molecu- 
lar number conservation. Furthermore the force balance 
condition. Equation ©, becomes 



2i?ocro,0 +Po - Po = 



(24) 
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which is Laplace's law with being the average pres- 
sures in the bulk. The value of the tension ctq ^ is in prin- 
ciple related to the value of a, but we will choose gn to 
be the free parameter here. The reader is referred to [ij 
for further details on the subject of membrane tensions 
and their conjugate areas. 

Calculating to first order in the deviations from a 
spherical shape is more tedious. Especially we need 
to find expressions for the hydrodynamic forces in 
terms of the state variables of the membrane. Fortu- 
nately, a solution to the Navier-Stokes equation at low 
Reynolds number in spherical geom etry is known. This 
solution is called the Lamb solution [ij, and we can use 
it to find in terms of Having worked this out 
we can then eliminate a^/m and from the dynam- 
ics of uim and nA.im to arrive at a closed set of first 
order differential equations for Uim and riA.im- n-^.^m 
will not enter the equations for ugm and nA,im, because 
9^y/(9nA5n5:)|„i=o = since the system should then 
be symmetric with respect to reversing the two sides of 
the membrane. We do not write out the intermediate re- 
sults of the full calculation here. Instead we refer to ^3 
where an analogous calculation was performed in detail. 
We state the final result by collecting the dynamical fields 
into a column vector lo^m, where the transpose of w^m 
is wj^^ = {uim,nA,im), and then write the differential 
equations as 



-Aw 



im 



(25) 



Here the matrix B represents the dissipative forces 



B 



?3 ^r+u'-i 




(26) 



and the matrix A the elastic ones 



A 



_A|o(^ + 2)(^-l) 



A a 

2 



a(^-|-2)(^-l) 





xRl 



where 



El = {K£{e + 1) + ao,^Rl) {e + 2){£-l) 



(27) 

(28) 
(29) 



no,s ,TiA=0 



III. THERMAL NOISE 



Equation (|25|l represents the near-equilibrium dynam- 
ics of the average shape of the membrane. To include 
the effect of the thermal forces on the d yna mics of the 
shape we follow the Langevin approach |l3l | and add a 
white-noise term to the dynamic equations. The result 
is the stochastic equation 



Win 



-B-^Awf 



'^thermal 



(30) 



where the statistics of the noise obeys 

(^thco„al(O^Lrmal(^')) = TthcrmaKJ {t - t' ) , (31) 

with T7jj^j,j.jjj^[ representing the complex conjugate of 
^thermal- Correlation matrix Fthermai can be deter- 

mined from the condition that the long time limit of 
the fluctuations calculated from Equation (|30l) should 
equal the fluctuations that can be calculated directly 
from the free energy using equilibrium statistical mechan- 
ics. The equilibrium statistical mechanical calculation is 
performed in appendix IXI with the result 



oq 



(32) 



To find the long time limit of the fiuctuations from Equa- 
tion (|30() we first note that Equation H30() has the formal 
solution 



(33) 

Squaring this and taking the average we find that the 
value at long times of the fiuctuations is 



WimW 



dt e-(^"'^)*r. 



thermal^ 



(34) 

Using B-^B^^A)"^ = {B-^A)B-^ it can easily be 
checked that the choice of rthermai that ensures the agree- 
ment between the two ways of calculating the fiuctuations 
is 



r 



thermal 



2kBTB-^ 



(35) 



IV. ACTIVITY 



Having set up the stochastic equations for the near- 
equilibrium dynamics we will now modify these to include 
the effect of the active force-dipoles. Most of the work 
needed to do this was done in 0, where it was found 
that the important contribution of the force-dipoles to 
the dynamics was to modify the force balance condition 
for the membrane to 



frs + /a 



0, 



(36) 



where f^^^ is the force per area generated by the activity. 
In this force was obtained as a divergence of a stress 
tensor 



/a 



where 



T. 



a(3 



a/3 



Hg 



af3 



-K 



Q 



(37) 

(38) 
(39) 
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The quantity adip is a modification of the tension induced 
by the activity. Due to the fast balancing of forces in 
the tangential directions of the membrane, fluctuations 
in (Tdip will immediately be counterbalanced by changes 
in tr^, and therefore only the average value ctq, dip will have 
an influence on the dynamics of the membrane shape. A 
discussion of the magnitude of this contribution can be 
found in 7]. Q is a measure of how much the bending 
moments in the membrane is modified by the activity of 
the force-dipoles. The model we will use for Q is 



Q = F{2)nA + S 



(40) 



The constant F(^2) is related to the parameters of Q by 
being the quadrupole moment of the shifted force-dipole 
of an active protein: 



F{2) - 



((■ 



(41) 



Here Fa is a force and w'^ and are length scales re- 
lated to how the active proteins insert themselves in the 
membrane. S represents the noise arising when the force- 
dipole generated by an individual protein is not constant 
during the cycle of the activity. It has zero average 



(S) = 



(42) 



and we model its autocorrelation function similarly to 
Q, 0, 0| by combining the self-diffusion of the pro- 
teins with an exponential decay of the correlations in the 
proteins internal state 

rano,sG,p(0, 0, 6', \t - i'|)e-l*-*'l/^-^ . (43) 

where Y a is a phenomenological constant, Tq a correlation 
time and Ggp is the single particle density correlation 
function for the active proteins. Ggp is the solution of 
the diffusion equation on a spherical vesicle 

^G,p(^, </), 0', 0', i) = i?,pi?„9"G,p(0, 0, 0', 0', t) , (44) 

with the initial condition 

G,p(^, 0', 0', t\^^ = -^5(6 - 9')5i<ly - 0') , (45) 
V9q 



where Dsp is the single particle diffusion constant and 
the spatial derivatives act on (9, (p). 

Note that, contrary to [EH) we have for simplicity 
ignored any dependence of the protein activity on the 
local curvature of the membrane, but included internal 
noise in the model instead. 

We calculate the stochastic equations of motion for utm 
and riA/m as before using the force balance condition 
of Equation (|36|l . instead of Equation jnj. This gives a 
changed set of stochastic equations which becomes 



-Cwp 



^thermal ' ''7a 



(46) 



with a modified dynamical matrix 

C = B-^{A + Aa) , 
where 

A _( ^7o.dipi?^(£-f2)(^-l) Ml(£-f2)(£-l) 
^" ^ ' 



(47) 



and an additional noise term 

Ro{e + 2)ie-l)Sem/2 



Va = 



where (Sem) — and 

{S*f>^{t)Sl'm'{t')) = 







(48) 



(49) 



l*-*'l/"^fe'<5w . (50) 



The time scale ti is a combination of and the relax- 
ation time for the single protein diffusion 



1_ _ 1 ^ Dspi{t+l) 



(51) 



V. FLUCTUATION SPECTRUM 

To find the fluctuation spectrum from Equation H46(l 
we first note that formula (|34|l is now modified to 



(52) 



where 



i — t thermal H ttt: I 1 — — I X 



1 




Rf, V 277i?3(4^3-f 6£2-l) 

-1 



C^ + -I 



C + -I 



I A V 1 







(53) 



If we assume that G has two distinct eigenvalues Ai and 
A2 then we can use the Caley-Hamilton theorem to find 
that we can write the exponential of —Ct as 



-Ct 



1 



-{(Aie-^=*-A2e-^^*)/ 



(54) 



A1-A2 

-f (e-^i* - e-^^*) G} 
Using this formula together with 

Ai + A2 = TrG, AiA2=dctG, (55) 
we can perform the integral in Equation 1)52(1 to find 



--"^-> -mZ^G (T^^ ^ - ^) r (TrG / - Cf 



1 



2TrG 



(56) 
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Equation H56|l is not singular at Ai = A2 and thus work 
even if the eigenvalues are identical. 

Using Equation H56(l we finally find the spectrum of 
the shape fluctuations including the contributions from 
the active force-dipoles 



{\uim\'^) 



AEi{l + T^i/TD,i) J ' 



(57) 



Here we have introduced the activity modified energy 
scales 



Ei^{Kl{l+l) + dRl){i + 2){l- 

z e 



1) 



(58) 

E^ = I^Ke{£ + 1) + a^RlJ (£ + 2){£ - 1) , (59) 

where the modified tension and bending parameters are 

(60) 
(61) 



~ e 
K 



<7 = G 



t + fo,dip , 

A(A-f(2)) 

4x 

A(A-J^(2)) 



(62) 



We have also introduced timescales for shape (r^.^) and 
protein density {td,i) relaxations 



Tk.1 



TD,i = 



■nRlju^ + 6£^ - I) 
EAe + 1) 



DAA£i£ + 1) 



where 



Daa = f^AAX 



(63) 
(64) 

(65) 



is the diffusion constant for nA- Finally we have intro- 
duced the fraction 



E + + Ei + ^Ee 



('^ 


^i^)e,^ 


-^EH 













(66) 



which determines the importance of the internal active 
noise. 

A plot showing how the force-dipolc activity modifies 
the fluctuation spectrum is shown in figure ^ The pa- 
rameters have been chosen in accordance with the esti- 
mates and findings of for a lipid membrane with the 
proton pump bacteriorhodopsin incorporated. Thus we 
have taken k = lOfceT with k^T ^ 4 • lO^^^ J. The 
protein potential V is taken to be purely entropic and 
thus X = k-QT/riQ^Y: with no.E = 10^® m^^ for the total 
density of proteins. We have assumed that the single 
protein and collective diffusion constants are identical 



. 0001 ^^i^rj. 

0.00001 - 

10-' : 

10"° 



5 10 20 50 100 

FIG. 1: The fluctuation spectrum, Equation 1571 . The pa- 
rameters are as explained in the text for the dotted curve. 
The solid curve is the equilibrium result, i.e. -^(2), Fa and 
o"o,dip has been set to zero. For the dashed curve we have 
taken Fa = such that only the internal active noise of the 
proteins is ignored. 
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FIG. 2; The fractions xi (solid curve), Tk,i/ti (dash-dotted 
curve), Tk,xI'^d,i (dashed curve) and Ti/To,e (dotted curve) 
as functions of I. 



and taken them to be -Daa = Ds-p — lO^^^m^/s. For 
simplicity we have ignored any asymmetry in the non- 
active interactions between the proteins and the lipids 
and taken A = 0. We have taken the vesicle to have a 
radius of i?o — 10 /im with a tension cro,0 — 5 • 10^* N/m 
which is increased by 50% when the proteins are active 
to (7 = 7.5-10"^ N/m. The relation between the constant 



F(2) here and the parameters of [6| is that F^ 



(2) = 2wVa, 

where w is the membrane thickness and Va has the 
units of a force-dipole. Using the values of this gives 
^(2) = 10 nm • K. The amplitude of the internal noise Fq 
should be related to f (2) Etnd here we have simply taken 
Va — F^2) ■ ^'-'^ characteristic time of the active pro- 
cess we have taken the duration of a bacteriorhodopsin 
photocycle, Tq = 5 ms. Finally the surrounding medium 
has the viscosity of water 77 — 10~^ Pa • s. 

Different ways in which the activity of the force-dipoles 
modify the fluctuation spectrum can be seen in flgure^ 
First of all the increase in the tension, tJdip, decreases the 
fluctuations at large wavelengths (small £). Secondly, the 
thermal fluctuations of riA enhances the shape fluctua- 
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tions at short wavelengths through the active contribu- 
tion to the bending moment: -F(2)?^a- This was the ex- 
planation given in 0| for the increase in the fluctuations 
observed in However, it can also be seen in the figure 
that the internal active noise can give a third modifica- 
tion at short wavelength, which sets in gradually as xi 
increases to become of order one, see figure El Note that 
eventually as the modenumber i goes to infinity, the con- 
tributions of the activity to the fluctuations scales as 
i.e. in the same way as the thermal fluctuations. In the 
classification scheme of 0, 0| the wavelength behavior 
of the fluctuations in the force-dipole model corresponds 
to the case of "curvature-force" . 

The physics behind the gradual onset of the effect of 
the internal active noise can be understood from the 
time-scales of the problem. Using the explicit values of 
the parameters chosen for the plot we find that at high 
£ we have Tk,j ~ 100s/£'^ and tdx — 100s/£^. Looking at 
the fractions of time scales in Equation H66I) we see that 
all of these fractions will be at maximum of the order 
unity except Tk,i/ti (see also figure|21l, which will be big- 
ger than one for i small enough such that r^/ is larger 
than Ta- Since this fraction enters in the denominator of 
Xf,^ we will have that xi is small for ^'s that are smaller 
than, say, 20. The physical origin of this is that when the 
membrane shape dynamics is slower than the dynamics 
of the internal active noise, then the shape cannot move 
fast enough to accommodate the temporary stresses that 
the internal active noise induces. Thus the noise will not 
influence the movement of the shape at long wavelength. 

Note that the timescale for diffusion to,i is slower than 
the timescale for membrane shape relaxation r„_£. A sim- 
plifying approximation in Equations l|56p and (|66() would 
thus be to discard terms proportional to Ti^j/tb/- 



VI. MICROPIPETTE EXPERIMENTS 

To understand what impact including the non-thermal 
noise represented by Ta in the force-dipole model has on 
interpreting existing experimental data, we will brieffy 
discuss micropipette experiments in this section [l9| . 

In micropipette experiments a relation between the 
area of a vesicle and the tension in its membrane is mea- 
sured. The vesicle is aspirated to a pipette at low suction 
pressure inducing a tension in the membrane of the size 
/j,N/m and above The tension is then increased by 
increasing the suction pressure and it is measured how 
the visible area of the membrane changes in response. 
At low tension it is usually assumed that the increase in 
visible area is taken from the excess area stored in the 
fluctuations of the membrane without the true area being 
increased. 

We can calculate the excess area, which we will denote 
by a' to distinguish it from the a of Equation ((Sj, for the 
force-dipole model applied to an spherical vesicle from 
Equation (|57|l . Expanding the area of the vesicle surface, 
A, to second order in u^m one finds for a spherical vesicle 



with fixed volume 
A 



- 1 



I- {t + 2){l-l)\ut. 



Stt 



|2 



(67) 



1=1 



where /max is a cut-off wavenumber beyond which the 
free energy of Equation is not expected to be valid 
anymore. Due to the relatively high tensions in the mi- 
cropipette experiments (in comparison with what can be 
achieved in video microscopy experiments) the fluctua- 
tions of the low wavenumbers H do not contribute much 
to the excess area. We can therefore simply replace fac- 
tors of, say, £ -I- 2 in Equations H57|) and (|67|l with / and 
convert the sum to an integral to simplify the calcula- 
tion. As mentioned at the end of the last section we 
can also put factors of Ti^^i/td/ to zero. To simplify the 
expression for xi we will simply take A to be zero such 
that Xf = l for high I, and assume that the tension is 
high enough such that the terms in the square brackets of 
Equation H57|) is suppressed by the tension for the lower 
I where xi is not unity. With these simplifications we get 



a 



log 



Sttk 



~aRl 



log 



"''max 



where 



= 1 



(68) 



(69) 



is an effective temperature whose ratio with T is a mea- 
sure of how much the excess area has been increased 
by the activity induced fluctuations in comparison with 
equilibrium. 

In ^, ^] a factor two to three increase of T^^ in mi- 
cropipette experiments on active membranes were re- 
ported. This increase was assigned to the averaged force- 
dipole activity represented by ^(2) Equation l|69|) . But 
as Equation H69|l shows, the noise of the activity repre- 
sented by Ta might well give a contribution of compara- 
ble size. Actually, with the estimated values chosen in 
the previous section the contributions of the correspond- 
ing two terms in Equation (|69|) become exactly equal. 
Further experiments are needed to make it possible to 
distinguish the contributions of the two mechanisms. An 
obvious candidate for such experiments is a video mi- 
croscopy experiment, where the two contributions might 
be distinguishable due to the onset at high ^-modes of 
the internal noise represented by Ta discussed in the last 
section. However, a limit of the video microscopy experi- 
ments is the resolution of the experimental setup, setting 
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a limit for the range of £-modes that can be obtained in 
such an experiment. Currently this limit seems to be the 
first 20 to 30 modes 0| , making the possibility of observ- 
ing an onset very sensitive to the actual parameters, k, 
Ta etc., of the system. 



VII. CONCLUSION 

In this paper the Gibbs formalism of 0, was used to 
calculate the quasi-spherical shape fluctuation spectrum, 
given in Equation lf77|) . for the force-dipole model of 
This calculation allows for further testing of the force- 
dipole through video microscopy experiments, something 
which is needed to test the model more stringently. The 
shortcoming of the micropipette experiments is that they 
only provide one data point for each set of experimental 
conditions. A video microscopy experiment on the other 
would provide a data point for each modcnumbcr i in the 
spectrum that can be measured. 

Following a general viewpoint promoted in an addi- 
tional contribution due to temporal noise in the strength 
of the force-dipoles of the individual active proteins was 
included in the calculation. The contribution of this noise 
was found to be insignificant at long wavelengths, but it is 
turned gradually on as the mode number i increases when 
the corresponding characteristic time scale for membrane 
shape relaxation crosses the time scale associated with 
the correlations in the noise from the activity. This grad- 
ual onset offers a possible way to distinguish the contri- 
bution from those of the original force-dipole model of 
when the shape fluctuation spectrum obtained by for 
instance a video microscopy experiment is analyzed. 
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APPENDIX A: EQUILIBRIUM FLUCTUATIONS 

For a column vector wgm with components that per- 
form small fluctuations about an average value zero we 
can calculate the equilibrium fluctuations as 



= kBTA- 



where the matrix A is the second derivative of the free 
energy 



A 



(A2) 



However, for the membrane case studied here there is a 
subtlety with respect to which fields to use in ivim- When 
we take the derivative of the free energy F in Equation 
Q we fix the fields y/griA and not . This is the reason 
that we choose for Wim (see also 10] ) 



(R-Ro) 



(A3) 



The transformation to the wgm we have used in the main 
part of the paper is then at linear order 



where 



L = 



I = Lwim , 

1 
1 

-2no,s 1 



(A4) 



(A5) 



The free energy. Equation Q, without the Lagrange 
multiplier term for but with a tension parameter CTq 
included instead to control the area is 

F = I dA [2KH^+AnAH + V{nA,ns)+(j'o] ■ (A6) 
Jm 

When we expand the free energy to second order in 
to find A we have to take care of the constraints on 
volume and number of molecules in the membrane. The 
volume constraint tells us that 



dA R-n 



(A7) 



M 



should be constant. To second order in we have 



v = m 



E 



(A8) 



This equation will be used to eliminate Mqo and express 
the free energy F in terms of the remaining u^ra when 
F is expanded to second order in Ui„i- Note that uoo 
depends on uim, ^ 7^ through terms which are second 
order in uim- We will therefore regard uqo as a second 
order term when we expand the free energy. 

The constraints on the number of particles, i.e. 
Jjyj dA UA is the total number of molecules of type A, 
simply constrains the last two components of tDoo to be 
zero. 

Using the constraints described above we can expand 
the free energy to second order in wim to find 



1 



F = constant H — 



(Al) where 



A 



Aflo(£+2)(£- 



-1) 



2 ^ 



2 

?2 



(A9) 



-2i?2x^no,s 



xR'o 







z?2 
X Ro 



Wirn—0 



(AlO) 
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and the tension that replaces (To,dip in the expression for 
El, Equation |(2HI), is 



where 



(7o = (j'q + V - 



'X^ is simply 



dV 



dri' 



dV 
duA 



"o,A=0,no,i: 



9n|, 



(All) 



(A12) 



ARo{e+2){e-i) 

2 





Aflo(^+2)(£-l) 
2 

?2 



xRo 





no,A=0,no,z; 



Using the linear transformation, Equation l|A4|l . we 
find that the fluctuations of wim with £ are 






p2 



(A14) 



WlmW 



(A13) This is the result used in Equation (|32|l . 
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